digamma(9) - 2.140641477955609996536345
digamma(2.5) - 0.7031566406452431872257
digamma(0.1) - (-10.42375494041107679516822)
digamma(7e-4) - (-1429.147493371120205005198)
digamma(7e-5) - (-14286.29138623969227538398)
digamma(7e-6) - (-142857.7200612932791081972)
digamma(2e-6) - (-500000.5772123750382073831)
digamma(1e-6) - (-1000000.577214019968668068)
digamma(7e-7) - (-1428572.005785942019703646)
digamma(-0.5) - (.03648997397857652055902367)
digamma(-1.1) - 10.15416395914385769902271
digamma(-1)
if digamma(0) ~= -Inf
  error('digamma(0) should be -Inf');
end
if ~isnan(digamma(-Inf))
  error('digamma(-Inf) should be NaN');
end
